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Abstract 

Context. Galaxy metallicities have been measured to redshift z~2 by gas-phase oxygen abundances of the interstellar medium using 
the R 2 3 and N2 methods. Galaxy stellar metallicities provide crucial data for chemical evolution models but have not been assessed 
reliably much outside the local Universe. 

Aims. We determine the iron-abundance, stellar metallicity of star-forming galaxies at redshift z~2, homogeneously-selected and 
observed as part of the Galaxy Mass Assembly ultra-deep Spectroscopic Survey (GMASS). 

Methods. We compute the equivalent width (EW) of a rest-frame mid-ultraviolet (mid-UV), photospheric absorption-line index, the 
1978 A index, found to vary monotonically with stellar metallicity by Rix, Pettini and collaborators (R04), in model star-forming 
galaxy (SFG) spectra created using the theoretical massive star models of Pauldrach and coworkers, and the evolutionary population 
synthesis code Starburst99. The 1978 A index is sensitive to Fe III transitions and measures the iron-abundance, stellar metallicity. 
To accurately determine the 1978 A index EW, we normalise and combine 75 SFG spectra from the GMASS survey to produce a 
spectrum corresponding to a total integration time 1652.5 hours (and a signal-to-noise ratio ~100 for our 1.5 A binning) of FORS2 
spectroscopic observations at the Very Large Telescope. 

Results. We measure a iron-abundance, stellar metallicity of log(Z/Z Q ) = -0.574 + 0.159 for our spectrum representative of a 
galaxy of stellar mass 9.4xl0 9 M Q assuming a Chabrier initial mass function (IMF). We find that the R04 model SFG spectrum for 
log(Z/Z Q ) = -0.699 solar metallicity provides the best description of our GMASS coadded spectrum. For similar galaxy stellar mass, 
our stellar metallicity is ~ 0.25 dex lower than the oxygen-abundance, gas-phase metallicity quantified by Erb and collaborators (E06) 
for UV-selected star-forming galaxies at z = 2. 

Conclusions. We measure the iron-abundance, stellar metallicity of star-forming galaxies at redshift z~2 by analysing the 1978 A 
index in a spectrum created by combining 75 galaxy spectra from the GMASS survey. We find that our measurement is ~0.25 dex 
lower than the oxygen-abundance gas-phase metallicity at similar values of galaxy stellar mass. We conclude that we are witnessing 
the establishment of a light-element overabundance in galaxies as they are being formed at redshift z~2. Our measurements are indeed 
reminiscent of the o--element enhancement seen in the likely progenitors of these starburst galaxies at low-redshift, i.e. galactic bulges 
and early-type galaxies. 

Key words. Methods: observational - Galaxies: abundances - Galaxies: evolution - Galaxies: formation - Galaxies: high-redshift - 
Galaxies: starburst 



1. Introduction 

Galaxies are believed to form by the merging of dark matter 
haloes, the isolated collapse of a gas cloud and subsequent gas 
accretion, or a combination of these scenarios. Photometric and 
spectroscopic surveys measure the rate of galaxy mass assem- 
bly and the variation in galaxy star formation rates with redshift 
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* Based on observations obtained at the ESO Very Large Telescope 
(VLT) as part of the Large Programme 173.A-0687 (the Galaxy Mass 
Assembly ultradeep Spectroscopic Survey). 



dLillv et all \]99a iMadau et all Il996t iGiavalisco et ail l2004at 
Cimatti et al. 200l- 



Galaxies are unlikely to evolve as closed boxes but ex- 
change gas, metals and star s with their environment. Following 
the early work of iLarsonl (1 19741) cosmological models pro- 
pose that mass outflows and winds driven by supernovae and 
AGN activity suppress star formation (via so-called "feed- 
back") increasingly effectively in galaxies of increasingly lower 
mass, and enrich the interstellar medium (ISM) with prod- 
ucts of stellar nucleosynthesis (|Springel & Hernquistl [2003; 
iBenson et all 120031: iDe Lucia et all 12004 ICroton et al] [2006; 
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iBower et al.ll2006t iBertone et al J 120071: iFinlator & Dave1l2007l) . 
Such winds are required for example to prevent the over- 
production of galactic low-mass stars and describe the ob- 
served metal-enrichment history of the Universe. In the local 
Universe and to a redshift of z~2 a "mass-metallicity rela- 
tion' 
ies 

2004; Savaglio et al 



is observed such th at more massive star - forming galax- 
ar e more metal-rich dLg queux et al. 1979; Tremonti et al. 

" 120051: lErb et all 120061: iLee et all 12006). 



Tremonti et al] d2004l) (hereinafter T04) found a tight correla- 



tion between the gas-phase metallicity and stellar mass of 53,400 
star-forming galaxies of redshift 0.005<z<0.25, using imag- 
ing and spectroscopic data from the Sloan Digital Sky Survey 
(SDSS). Their relation has a lcr scatter from the median of 
only 0.1 dex and is linea r for stellar masses 8.5 < log(M/M ) 
< 10.5. At higher redshift IFinlator & Pavel d2007l) were able to 
more accurately reproduce the gas phase mass- metallicity re- 
lation at z~2 using "momentum-d riven winds" (iMurrav et alj 
2005; Oppenh eimer & D ave 2006) instead of constant veloc- 
ity winds traditionally proposed. "Momentum-driven winds" re- 
move metals into the intergalactic medium by radiation pressure 
as photons are absorbed and scattered by dust, and supernovae 
and galactic winds. Unlike thermal energy, momentum cannot 
be radiated away and "momentum-driven winds" can transport 
metals to significant distances outside the galaxy. 

Measuring galaxy metallicity as a function of redshift 
probes the rate of chemical enrichment. In the local Universe 
metallicity is meas ured using both spec troscopic absorption 
and emission lines. Gallaz zi et al.l d2005l) found good agree- 
ment between stellar and gas-phase metallicities using SDSS 
galaxy spectra for which both absorption and emission lines 
were measurable. Faber and collaborators at Lick Observatory 
pioneered measurement of luminosity-weighted stellar age 
and metallicity of nearby early-type galaxies by develop- 
ing the "Lick system" of spectroscopic absorption-line in- 



dices (Burstein et al. 


1984t iFaber et alj 119851: iGonzalez 119931: 


Gorgas et alJ 11993: Worthev et al.l 


1994; Worthev & Ottaviani 


1997; iTrager et alJll998l) for the rest-frame optical wavelength 



range ~4000-6200A, a nd the single-bu rst evolutionary popula- 
tion synthesis models of Worthev ( 1994). The "Lick system" was 
further developed b vlTrager et al.l (1 2000) and Thomas, Maraston 
and collaborators dThomas et al.l l2.Q03 . 2004) to account for the 
effects of light-element overabundance ratios on Lick/IDS in- 
dex measureme nts using respectively the empirical r elations of 
iTripicco & Belli d 19951) . a nd the theoretical models of lMarastonl 
d2005l) and calibrations of iKorn et all d2005h . 

At high-redshift emission lines are more readily de- 
tected than absorption lines and galaxy metallicity has been 
measured using primarily the oxygen-abundance, gas-phase 
metal licity of the IS M of star- forming galaxies, using the 



dPagel et alJll979h and N2 dStorchi-Bergma nn et all 
Rai mann et alj l2000j) me th ods dKobulnickv & Zaritsli I 
Kobulnickv & Kod l2000l: ICarollo & Lilly! 12001b iLillvi 
20031: iKobulnickv & Kewleyl 12004 IShaplev etail [2004; 
Maier et all 120051: lErb et all I2006T lLamareille et alj [2006; 



et al. 



Maie r et al.ll2006l) . At z>l the gas m etallicity of the ISM has 



been constrained using the CRM 1550 (Mehlert et al. 2002) and 
Mgll and Fell absorption lines dSavaglio et alj 120041) . Pettini 
and collaborators have measured the chemical abundances of 
the ISM using rest-frame UV absorption lines in a spectrum 
of the gravitational ly-lensed Lyman-g galaxy MS 1512-cB58 
at redshi ft z~2.73 ( Pettini et al. 2002b) and damped Lyman-a' 
systems (Petti ni et alj l2002ah . Stellar metallicit y data are re- 
quired to fully test chemical evolution models, de Mell o et all 
(2004) provided a first constraint of stellar metallicity for a 



small sample of very massive (> IO^'Mq) z ~ 2 galaxie s from 
the K20 survey dCimatti et alj|2002h iDaddi et al.ll2004al) . They 
showed tentative evidence for solar or supersolar metallicity but 
with l arge unce r taintie s due to limitations in the achieved S/N 
ratio. iRix et al] d2004l) compared the iron-abundance, stellar 
metallicity for MS 1512-cB58 at z=2.73, and Q1307-BM1163 
at z= 1.411, with previously determined light-element abun- 
dances (magnesium, silicon and sulphur) and iron-abundance 
gas-phase metallicities of the ISM. For MS 1512-cB58, the 
iron-abundance, stellar metallicity is a factor of two higher 
than the light-element abundance gas-phase metallicities. The 
iron-abundance gas-phase metallicity is a factor of four lower 
than the light-element abundance metallicities: this is attributed 
to dust depletion of the ISM and the time delay in release of 
products of supernovae la (Pett ini et alj l2002~bb . For Q1307- 
BM1163, the iron-abundance, stellar metallicity is a factor of 
two higher than the oxygen-abun dance, gas - phase metallicity 
measured using the N2 method. lErb et al] d2006l) compared 
two rest-frame mid-UV spectra (obtained for four of the six 
bins of stellar mass examined in their mass-metallicity relation 
analysis i.e. two lowest and two highest of six bins of mass) with 
star-forming galaxy s pectra g enerated u s ing th e Starburst99 
dLeitherer et alj 1 19991) model. lErb et alj d2006l) found good 
agreement between the gas-phase metallicity oxygen abundance 
determined using the N2 method and the metallicity of the 
model spectrum that best described the rest-frame mid-UV 
spectra as judged by eye. 

We present the first robust measurement of iron-abundance, 
stellar metallicity for mid-IR selected star-forming galaxies 
(SFGs) at z~2. We measure a ph otospheric absorption-line in- 
dex, defined by IRix et all d2004l) (hereinafter R04) in the rest- 
frame mid-ultraviolet galaxy spectrum (mid-UV) that is domi- 
nated by the light of massive stars. R04 modelled star-forming 
galaxy spectra usin g the evolutio n ary p opulation synthesis mod- 
els Starburst99 dLeitherer et alj d 19991) : hereinafter L99) and 
non-local ther mal equilibrium (NLTE ) stellar atmosphere mod- 
els of OB stars ( Pauldrach et alJl200ll) . Using stellar atmosphere 
models, the wavelength coverage of theoretical star-forming 
galaxy (SFG) models was extended redward to ~2100 A. Star- 
forming galaxy model spectra were generated for five values of 
metallicity, 0.05, 0.2, 0.4, 1.0, and 2.0 solar, assuming a contin- 
uous star formation history and a Salpeter IMF. R04 found that 
the equivalent width (EW) of an absorption-line system, centred 
on 1978 A within an Fe III transition blend, varied monotoni- 
cally with metallicity, after 100 Myr of star formation. They de- 
fined the EW of the "1978 A index", EW(1978), to have limits of 
1935 A and 2020 A avoiding the Al III /U1855, 1863 interstel- 
lar lines, the nebular emission line C III]/U909 and other weaker 
lines. 

We measure the 1978 A index in a spectrum created by 
adding 75 star-forming galaxy spectra from the Galaxy Mass 
Assembly ultradeep Spectroscopic Survey (GMASS) data corre- 
sponding to over 1652.5 hours of integration time with FORS2 at 
the VLT. Given the difficulty of measurement of the 1978 A in- 
dex and the limited range of galaxy stellar masses, we do not 
analyse spectra as a function of stellar mass. It would indeed 
be challenging to achieve this measurement for massive galaxies 
alone, because the m ost massive galaxies are fainter in the rest- 
frame mid-UV (e.g. IDaddi etai1l2004al: Hong et al.ll2006l) and 
they are of course rarer (cf. Figure|2|. 

In Section [2] we provide a brief description of the GMASS 
survey. In Sections [3] and [4] we describe respectively, the 
selection of, and the combination of GMASS star-forming 
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galaxy spectra. In Section [5] we outline the measurement of 
the 1978 A index and its error, and contrast our observed data 
with the R04 star-forming galaxy model spectra. In Section [6] 
the 1978 A index measurement is compared with predictions of 
galaxy stellar and gas metallic ity from the cosmological simula- 
tions of lFinlator & Pavel (120071), and t he oxygen-abundance gas- 
phase metallicities of lErb et all d2006l) and we briefly discuss our 
results and conclusions. 

We use AB photometric magnitudes, and assume a WMAP 
cosmology with Q a,Qm - 0.73,0.2 7, and h =Ho[km s _I 
Mpc^J/lOO = 0.71 (Spergel et alj2003l). W e assume a Chabrier 
initial mass function (IMF) dChabrierl2003l) from 0.1 to 100 M B . 
Whether the IMF is universal between galaxies and at different 
redshift is a subject of ongoing debate. The empirical measure- 
ment of the Galactic IMF is not well determined and the uncer- 
tainty increases for the individual ki nematical com ponents (disk, 
bulge and thick disk) of the Galaxy (Chabrier 2003), and to high- 
redshift. For a review of the observational constraints provided 
by measurement of the mass-to-light ratios of nearby galaxies, 
the redshift evolution of the F undamental Pla ne and galaxy clus- 
ter, element-abundances, see Renzini (2005). In this paper, the 
Chabrier IMF is chosen to provide a reference for all data shown 
and our conclusions are independent of this choice. 

2. The GMASS Survey 

GMASS ("Galaxy Mass Assembly ultra-deep Spectroscopic 
Survey"Q) is an ESO VLT Large Program project based on data 
acquired using the FORS2 spectrograph. A complete description 
of GMASS photometric observations, source detection, target 
selection, determination of photometric redshifts, mask prepara- 
tion and spectroscopic observations and data reduction are pro- 
vided by Kurk et al. (in preparation). 

The project's main science driver is to measure the physical 
properties of galaxies at redshifts 1.5 < z < 3 a critical range 
in the mass assembly of massive galaxies. Spectroscopy allows 
the measurement of reliable galaxy redshifts, stellar masses, star 
formation rates and metallicities. The uniqueness of GMASS is 
its 4.5jt/m selection that provides two major benefits: (1) it de- 
tects the peak of the stellar SEDs (for A resl = 1 .6//m) redshifted 
into the 4.5yum band for z > 1.4, and (2) it is sensitive to the 
rest-frame mid-IR emission, i.e. to stellar mass, up to z«3. The 
stellar mass completeness limits are log(M/M ) « 9.8, 10.1, and 
10.5 for z — 1.4, z = 2, and z = 3, respectively (for m.4.5 < 23.0 
and a Chabrier IMF). These mass limits allow mass assembly to 
be tracked for the precursors of today's massive galaxies to the 
most massive galaxies at z > 1.4. 

GMASS target selection was completed for the GOODS- 
South fielc0 from a region of 6.8 x 6.8 arcmin 2 chosen to match 
the FORS2 spectrograph field-of-view. All sources with Spitzer 
Space Telescope + IRAC data (Dickinson et al. in preparation), 
were selected to a limiting magnitude of m.45 < 23.0 (2.3 //Jy). 
Finally a cut in photometric redshift of z p hot > 1 -4 and two cuts 
in the optical magnitudes (B < 26.5, / < 26.5) were applied. 

Spectroscopic observations were completed for integration 
times of up to 32 hours, using either the blue 300V grism (for ob- 
served wavelengths 4000-6000 A), or the red 3001 grism (6000- 
10000 A) depending on target photometric properties. The spec- 
troscopic slit width was 1 arcsecond. Overall the GMASS sur- 
vey had an excellent spectroscopic redshift measurement suc- 
cess rate of ~85%. For objects without GMASS spectroscopic 

1 http://www.arcetri.astro.it/~cimatti/gmass/gmass.html 

2 http://www.stsci.edu/science/goods 



redshifts we adopted data from the literature where available, 
or photometric redshifts derived fro m our analysis using th e 
publicly available optical (HST+ACS) dGiavalisco et al.l20 04b), 
near-IR (VLT+ISAAC) (Retzlaff et al. in preparation) and IRAC 
images of the GOODS-South field (Dickinson et al. in prepara- 
tion). 

3. Selection of galaxy spectra 

Specifically for the present analysis, all spectroscopic data col- 
lected for the GMASS survey were visually classified as passive 
galaxy, star-forming galaxy or stellar spectra, or (if a spectrum 
had too low signal-to-noise) to have an uncertain spectral type. 

A spectrum was introduced into our analysis if classified 
as a star-forming galaxy (SFG) spectrum and if its rest-frame 
wavelength range included 1700-2100 A. In Figure|2]we present 
the distribution of redshifts, K$ apparent magnitudes and galaxy 
stellar masses of the 75 GMASS spectra selected. In Figure Q] 
we show a typical SFG spectrum created by simply average- 
combining all 75 GMASS SFG galaxy spectra chosen. 

In Table [1] we indicate the median physical properties of the 
75 selected galaxies i.e. for photometric data from the GOODS 
project, the apparent B, V, I and z photometric magnitudes de- 
rived using HST ACS imaging for the filters F435W (h goo ds), 
F606 W (v^,), F814W (w ,,), and F850LP (z goods ) respec- 
tively (Giavali sco et alj|2004bl) . the J, H, K$ photometric magni- 
tude, and absolute B, I and K magnitude; photometric and spec- 
troscopic redshift; the best-fit galaxy age in Gyr, star formation 
rate in solar masses per year, dust extinction in magnitudes, and 
stellar mass in solar masses. Galaxy stellar mass, age, star for- 
mation rate, and dust extinction were meas ured by fitting the 
evolut ionary population synthesis models of Bruzual & Chariot! 
(2003) to galaxy spectral energy distributions SEDs inferred 
from photometric data gathered for GMASS. The model as- 
sumptions were a Chabrier IMF, solar metallicity, and a star for- 
mation rate (SFR) proportional to e~ where r is the star for- 
mation characteristic timescale, and t the time elapsed since the 
onset of the current episode of star formation (Pozzetti et al. in 
preparation). 

4. Coaddition of star-forming galaxy spectra 

To be able to measure the 1978 A index and the galaxy iron- 
abundance, stellar metallicity, we carefully normalised and com- 
bined all 75 SFG spectra selected in Section[3] 

For all observed wavelengths, the flux-calibrated intensity of 
each spectrum was deredshifted to the galaxy rest-frame using 
the GMASS spectroscopic redshift, and then divided by the me- 
dian value of flux. A spectrum was created for the wavelength 
range 1700-2100 A by fitting a natural spline function at regu- 
lar intervals of 1.5 A. A linear function was then fitted to and 
divided into this spectrum to take out a remaining slope in the 
spectral shape. 

Due to the wide range of signal-to-noise ratios of the 75 SFG 
spectra (because of the diverse range of galaxy magnitudes and 
colours) and to ensure the most precise 1978 A index measure- 
ment, we optimised the signal-to-noise ratio of the coadded spec- 
trum using weights. Each spectrum was "weighted" by a mea- 
sure of its signal-to-noise ratio and, as a function of wavelength, 
by the relative contribution of background night sky to the spec- 
trum flux. To assess the weight due to signal-to-noise ratio of 
each normalised spectrum, we calculated the semi-interquartile 
range (s.i.q.r.) of flux in the wavelength range 1920-2050 A and 
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Table 1. The median galaxy physical properties of 75 GMASS spectra combined to produce a spectrum used to measure the 1978 
A index. Data presented are apparent B, V, I and z, J, H, Ks photometric magnitudes, absolute B, I and K magnitude (absB, absl and 
absK), photometric (z p hot) an d spectroscopic (z spec ) redshift, best-fit galaxy age in Gyr, dust extinction Ay, star formation rate in 
solar masses per year, and stellar mass in solar masses. 



b goods Vgoods hoods Zgoods J H Ks absB absl absK z phot z spec age(Gyr) A v SFR(M /yr) Mass(Mp) 

2463 2441 2402 23X2 2334 23T6 ZL87 ^2144 ^2L96 ^22~29 L91 L88 0T8 O80 26758 8.32x10" 




2xl0~ 19 ' 1 ^ J 1 1 L — ' k ■■ 1 1 L^-i i lJj l_j i_j i L_i i j, i i_u u i i I l_j i; i 

1300 1400 1500 1600 1700 1800 1900 2000 2100 

Rest — frame wavelength (A) 

Figure 1. We provide an average-combined spectrum of the 75 GMASS star-forming galaxy spectra selected for the present analy- 
sis. Dotted lines indicate interstellar rest-frame mid-UV absorption lines; dot-short dashed lines mark photospheric absorption lines; 
and the dot - long dashed line labels CIII]/i 1908.734 emission. Where two lines are close in wavelength and correspond to a single 
element, one label is shown. The R04 1978 A index equivalent width wavelength range is delineated by straight solid black lines at 
1935 A and 2020 A. 




spectroscopic redshift Ks SB apparent magnitude log 10 (M*) 



Figure 2. The distribution of GMASS spectroscopic redshifts, GOODS K s magnitudes and GMASS stellar masses of the 75 
galaxies for which GMASS star-forming galaxy spectra are coadded in this analysis. 
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then its quadrature, and subtracted in quadrature a term (0.03) es- 
timated to represent the contribution of intrinsic spectrum shape 
to the calculated s.i.q.r. value. A background sky spectrum was 
obtained by average-combining rows containing predominantly 
sky spectra in a CCD galaxy spectrum image. At all wavelengths 
the background sky flux was divided by its median value. The 
weight applied to each individual spectrum during combina- 
tion was the product of the inverse of this final background sky 
spectrum, and the inverse of the weight due to spectrum signal- 
to-noise ratio calculated as described above. A stellar mass of 
9.4xl0 9 M was calculated for the final coadded spectrum by 
weighting the GMASS stellar mass for each of the 75 galaxies 
by the identical weights applied to each spectrum during spec- 
trum combination. 



5. Stellar metallicity measurement and error 

The 1978 A index has a typical EW of 2-7 A (iRix et al.l2004 and 
its measurement is easily affected by a few deviant pixels in par- 
ticular inside the pseudo-continua wavelength intervals defined 
by R04. In this Section we describe the measurement of the stel- 
lar metallicity and its error, and assess the solidity of our results 
using independent methods. 

Firstly however we consider the appropriateness of the R04 
predictions for our observed data. The R04 empirical calibra- 
tion between 1978 A index EW and metallicity (Equation 8 of 
R04) is defined for specific model assumptions and a particular 
spectral resolution. The R04 models are applicable to galaxies 
that have been undergoing bursts of star-formation for > 100 
Myr, a reasonable assumption for star-forming galaxies at red - 
shift 2 dDaddi et al.ll2004bHShaplev et al.ll2004TlErb et al.ll2006l) . 
The rest-frame spectral resolution of the R04 model predictions 
is 2.5 A; this differs significantly from our rest-frame spectral 
resolution of 3.8 A FWHM. 

To be able to reliably measure stellar metallicity we derived a 
new empirical calibration between the 1978 A index and metal- 
licity for our rest-frame spectral resolution of ~3.8 A FWHM. 
To bring each R04 model spectrum to our resolution, we con- 
volved each spectrum with a Gaussian function of width 2.86 A 
FWHM. We measured the 1978 A index using an approach iden- 
tical to the method described for our observed spectrum below 
in Section IBTTI 

Our revision of the R04 calibration is identical in form to 
Equation 8 of R04 i.e. 

log(-) = C-£W(1978) + L> (1) 

Z 

but our values of the coefficients C and D are different i.e. for 
an EW(1978) < 6.0 A, we find C = 0.37 and D = -1.95. In 
Figure [3] we directly compare the original R04 1978 A index - 
metallicity calibration, and our revised calibration. This figure 
illustrates that for a given EW(1978) measurement, our revised 
calibration implies a higher value of metallicity than the original 
R04 calibration. We use our revised calibration to determine the 
galaxy iron-abundance, stellar metallicity. 

We remark that to ascertain our 1978 A index EW and its er- 
ror and complete Figures|H|5]and|6] a small positive offset of 1 A 
was applied to the dispersion axis of the coadded spectrum. This 
offset was necessary to optimise the match between our coad- 
ded spectrum and the R04 model spectra i.e. the wavelengths 
of the faint absorption-line spectrum inside 1935-2020A. The 
offset is similar to our wavelength calibration uncertainty and 




z/z e 

Figure 3. Our revised 1978 A index - metallicity calibra- 
tion (solid line) is compared with the original R04 relation 
(long dashed line). Our 1978 A index measurement and error 
(3.72+0.43 A) are indicated by a horizontal solid and two hor- 
izontal dotted lines respectively. Filled circle symbols indicate 
the 1978 A index measurements for each of the five R04 model 
spectra broadened to our rest-frame spectral resolution and their 
five values of stellar metallicity. Unfilled circle symbols show 
the corresponding metallicity implied by the original R04 rela- 
tion for the 1978 A index values measured for each broadened 
spectrum. 

corresponds to a sub fraction of a pixel in the GMASS spectro- 
scopic data. If the offset had not been applied our 1978 A index 
measurement would have changed by far less than the error we 
estimate. 

5.1. The 1978 A index measurement 

We measured the 1978 A absorption-line index equivalent width 
(EW) as defined by R04 in our coadded spectrum using fixed 
limits of 1935 A and 2020 A and by fitting the continuum using 
6 of the narrow (typically 3-5 A wide) pseudo-continua intervals 
for the wavelength range 1700-2100 A. For our 1978 A index 
EW measurement of 3.72 A our revised calibration implies an 
iron-abundance, stellar metallicity of 0.267 solar. 

5.2. Fitting R04 model galaxy spectra 

We independently determined the stellar metallicity by fitting the 
five R04 spectra (for metallicities of 0.05, 0.2, 0.4, 1.0, and 2.0 
solar), broadened to match our rest-frame spectral resolution, to 
our coadded observed spectrum for the wavelength range 1935- 
2020 A. In our so-called "renormalised method" we decided the 
best-fit R04 model spectrum by minimising the median of the 
quadrature of difference in flux between each broadened R04 
model spectrum and our coadded spectrum. We minimised the 
median flux difference by multiplying the model spectrum by 
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numbers between 0.9 and 1.1. We searched for the minimum 
of the quadrature of the median and not of the average because 
(i) the median is less sensitive than the average to deviant pix- 
els, and (ii) we could not reliably measure error as a function of 
wavelength in our coadded spectrum (although to first order this 
should be fairly constant for our wavelength range of interest 
here). In our direct method we fitted each broadened R04 model 
spectrum without any multiplication by a constant. In Figure 
[5] we plot the quadrature of the median difference in flux for 
the two different fitting methods at each R04 model metallicity. 
Using both methods we find that the Z = 0.2 solar R04 broadened 
model spectrum provides the best description of our coadded 
spectrum. This is in good agreement with the iron-abundance, 
stellar metallicity measured using the 1978 A index. 

5.3. Robustness of the measurement and error 

We assess the integrity of our stellar metallicity measurement 
and its error in different ways. 

On close examination of Figure [5] it is seen that the faint 
absorption lines of our coadded spectrum are well reproduced 
by the R04 model spectra. This provides reassurance that the 
model spectra, broadened to match our spectra resolution, can 
be used to measure the galaxy stellar metallicity of our data. 

The agreement between the two independent, best-fitting 
methods, completed in Section |5T2l implies that we are correctly 
and consistently assessing the continuum level of the coadded 
spectrum. This is critical because incorrect determination of the 
local (pseudo-)continua levels can be the main source of ran- 
dom and systematic error in the measurement of absorption-line 
equivalent widths. 

We further remark that when applying the renormalised 
method in Section 15.21 the constant required to minimise the 
flux difference between coadded spectrum and model spectra (in 
cases of the 0.05, 0.2 and 0.4 solar metallicity models), modified 
the continuum level of the model spectrum by typically 0.6%. 
The coadded spectrum continuum level is being reproduced to 
at least this level. A conservative 1978 A index error measure- 
ment would be ~0.5 A (the error corresponding to a shift of 0.6% 
in the level of the normalised coadded spectrum continuum). 

In Section 15.21 and Figure |4] the quadrature of the median 
difference between the coadded spectrum and its best-fit model 
spectrum, corresponds to ~1.2% of the coadded spectrum con- 
tinuum level. This implies that the assessment of the contin- 
uum in the coadded spectrum (normalised using 6 R04 pseudo- 
continua intervals) is accurate to a level of ~0.5%. Conversely 
the R04 model spectra are proven to reproduce an observed star- 
forming galaxy spectrum at the level of ^1%. 

Using the above results we conservatively assume that the 
noise of each bin of our coadded spectrum is ~1% (implying a 
signal-to-noise ratio S/N of 100 per 1.5 A rest-frame resolution 
element). 

We determined the 1978 A index error by simulating galaxy 
spectra affected by noise. As previously noted a major source 
of error in the 1978 A index EW measurement is inaccurate as- 
sessment of the continuum level using local pseudocontinua ei- 
ther side of the main index passband. 1000 noise-affected galaxy 
spectra were created by adding random noise to the mean flux 
in each of the 6 R04 pseudocontinua wavelength intervals. The 
1978 A index was measured in each spectrum as described in 
Section IO The standard deviation of the 1000 1978 A index 
measurements was 0.43 A. This was taken to be our index error 
and is very close to our independent assessments above. Using 
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Figure 4. For the wavelength range 1935-2020 A we plot the 
quadrature of the median of flux difference between our coadded 
spectrum and each R04 model spectrum broadened to match our 
rest-frame spectral resolution. These differences were calculated 
as described in Section l5T2l The median is shown as a function of 
R04 model spectrum metallicity in units of solar metallicity for 
both the direct and renormalised methods. For both methods the 
difference between our coadded spectrum and model spectra is 
a minimum for the R04 model spectrum of 0.2 solar metallicity. 
The dotted vertical lines indicate our stellar metallicity measure- 
ment and error derived in Section lBTTI 

error propagation in Equation [T] this error corresponds to an un- 
certainty in metallicity of 0.159 dex: this is the metallicity error- 
bar plotted for our data in Figure [6] 

Our measurement of stellar metallicity, derived using the 
1978 A index, and its corresponding error is 0.267 ^'g^ solar 
(log(Z/Z ) = -0.574 + 0.159). 

6. Discussion 

In Figure [6] we compare our measurement of galaxy iron- 
abundance, stellar metallicity, with the oxygen-abundance, gas- 
phase metallicity data of UV-selected star-forming galaxies from 
Erb et al. (2006, hereinafter E06), and the theoretical predic- 
tions of both galaxy gas-phase and stellar metallicity for dif- 
fe rent galaxy stella r mas ses from the cosmological simulations 
of Finlator "& Pavel d2007f) (hereinafter FD07).The esti mated er- 
ror of the N2 method itself as indicated in Figure 3 of Erb et al.l 
(2006) is approx imately 0.05 dex. The solar oxygen abundance 
measurement of Asplu ncTet al.l d2004l) has been subtracted in 
log-space from the oxygen-abundance gas-phase metallicities of 
E06; the stellar masses presented in E06 were calculated by as- 
suming a Chabrier IMF. 

Before venturing into interpreting Figure|6]we need to clarify 
that "metallicity" has a different meaning as measured by us, ob- 
served by E06 or calculated by FD07. In our case, the 1978 A in- 
dex measures a photospheric, Fe III absorption line of massive 
OB stars, i.e., of the most recent stars to have formed out of the 
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Figure 5. Plots of optimised best-fitting of our coadded spectrum with the five R04 model galaxy spectra broadened to match our 
spectral resolution. In each plot our coadded spectrum is shown by a black solid line. Broadened R04 model spectra are overplotted 
with red lines for metallicities = 0.05, 0.2, 0.4, 1.0, and 2.0. Black long-dashed vertical lines indicate the wavelength interval 
corresponding to the R04 1978 A index, and black solid horizontal lines delineate a normalised spectrum flux level of 1.0 . 



ISM. The E06 measurement instead refers to the gas-phase oxy- 
gen abundance for UV-selected star-forming galaxies using the 
N2 method. Finally, the FD07 "metallicity" refers to the mass- 
averaged abundance of all the stars formed during the whole pre- 
vious history of each galaxy, and corresponds to all heavy ele- 
ments having assumed a global yield (y = 0.02). Thus the three 



metallicities refers to three different galactic components, and to 
three different elements or combination of elements. Moreover, 
the three metallicity determinations are derived from completely 
different procedures, each affected by its own different system- 
atic errors. With this caveat in mind, we avoid overinterpreting 
the differences seen in Figure|6] 
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The ~0.2 dex offset between the E06 oxygen abundance of 
the ISM and the average metal abundance of the stars of the 
FD07 simulated galaxies is qualitatively in agreement with the 
expectation. Indeed, the ISM metallicity results from the whole 
past metal enrichment, while the average stellar metallicity in- 
cludes all previous, more metal poor stellar generations. 

More intriguing is the difference between our estimate and 
that of either the ISM (~0.25 dex) or the stars in the simulated 
galaxies (~0.05 dex). Since we measure the metallicity in mas- 
sive OB stars, one would expect it to be nearly identical to that 
of the ISM, given the short time elapsed since these stars formed 
(less than ~10 Myr). However, assuming that this difference is 
real and not the effect of systematic error in the different meth- 
ods, we must consider that the ISM metallicity is a measure of 
oxygen abundance, and our stellar metallicity is a measure of 
iron. As is well known, oxygen is produced only by massive 
stars exploding as Type II supernovae within at most few 10 Myr 
after formation, while iron is also produced by Type la super- 
novae that are characterized by an extremely wide distribution 
of formation-to- explosion delay times, from ~ 10 My r to more 
than 10 Gyr (e.g. lGreggidl2005[|Mannucci et alj|2005l) . 

This different time behaviour of the two supernova types is 
understood to be responsible for the a element enhancement (in- 
cluding the abundance of oxygen) seen in stellar systems that 
formed over a time short enough to correspond to a fraction 
of the SNIa distribution of delay times. Indeed, such a ele- 
ment enhancements are observed in the metal poor stars o f the 
Galactic halo (e.g. IWheeler et alj|1989t [Cavrel et all 120041). in 
metal poor and rich stars of the Galactic bulge (IZoccali et all 
2006; Fulbright et al. 2007), as well as in elliptical galaxies (e.g. 
Davies et al.ll993l:lKuntechnerll998tlHallidavll999l:lKuntscnnerj 
2000: iKuntschner et al.ll200lh iMorelli et alj|2004t iThomas et alj 
2005). Since the star-forming galaxies we are studying are un- 
dergoing a starburst phase, the likely progeny of many if not 
most of them are bulges and early-type galaxies which in the lo- 
cal universe exhibit an a enhancement, comparable to the 0.25 
dex offset seen in Figure[6] Thus, taking this figure at face value, 
it suggests that we may be witnessing the a enhancement being 
established among these galaxies. 

The original low-redshift mass-metallicity relation of T04 
lies ~0.5 dex above the relation of E06. E06 also remeasured the 
oxygen-abundance, gas-phase metallicity using the N2 method 
for the same 53,400 galaxies studied by T04. They found a sim- 
ilar offset although narrower at higher stellar mass close to solar 
metallicity as the [Nil] emission line used in the N2 method be- 
gins to saturate. 

We note that the weighting scheme implemented, and the re- 
quirement of sufficient signal in the rest-frame UV to be able to 
measure a galaxy redshift, is biasing our measurement against 
the most reddened z — 2 galaxies (and thus likely the most 
metal-rich). Although the long GMASS integration times of 15- 
20 hours are limiting this effect, still it might be expected that 
our result could be somewhat underestimating the true average 
metallicity of star-formi ng galaxies at t he probed stellar masses 
(we note that the data o flErb et all (2006) probably suffer an even 
stronger bias because of their UV-selection). 

7. Conclusions 

We have measured the iron-abundance, stellar metallicity of 
star-forming galaxies at redshift z~2, for a spectrum created 
by adding together 75 star-forming galaxy spectra from the 
GMASS survey. The iron-abundance, stellar metallicity is de- 
termined by measuring the EW of a rest-frame mid-UV, pho- 



tospheric absorption-line index, the 1978 A index, defined by 
Rixetal. (2004) using the theoretical massive star models of 
iPauldrach et alj d200ll) . and the evolutionary population syn- 
thesis code Starburst99. We measure an iron-abundance, stel- 
lar metallicity of 0.267 solar (log(Z/Z ) = -0.574 + 
0.159). This is lower by ~0.25 dex than the oxygen -abundance, 
gas-phase metallicity measured by lErb et al] d2006l) for similar 
galaxy stellar mass and redshift. At least part of this difference 
may be the result of different systemat ic errors betwee n the two 
estimates although both ourselves and Erb et al. (2006) find that 
systematic errors are very unlikely to reproduce our observed 
data results. We postulate that our measurements are reminiscent 
of the c-element enhancement seen in the likely progenitors of 
these starburst galaxies, i.e. galactic bulges and early-type galax- 
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Figure 6. The blue hexagon filled symbol indicates our measurements of galaxy stellar mass and metallicity. Red circle unfilled 
and black circle filled symbols pro vide cosmological model predictions of galaxy stellar and gas metal licity respectively at redshift 
z~2 from iFinlator & Pavel (120071) . Black triangle filled symbols are the gas-phase metallicity data of lErb et akl (120061) at redshift 
z~2. 
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